Coherence Window in the dynamics of Quantum Nanomagnets 
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Decoherence in many solid-state systems is anomalously high, frustrating efforts to make solid- 
state qubits. We show that in nanomagnetic insulators in large transverse fields, there can be a fairly 
narrow field region in which both phonon and nuclear spin-mediated decoherence are drastically 
reduced. As examples we calculate decoherence rates for the _Fe-8 nanomolecule, for Ni particles, 
and for Ho ions in LiHoxYi-zF4. The reduction in the decoherence, compared to low field rates, can 
exceed 6 orders of magnitude. The results also give limitations on the observability of macroscopic 
coherence effects in magnetic systems. 

PACS numbers: 



Decoherence has emerged as a major challenge, both 
for fundamental physics, and for attempts to make solid- 
state qubits (lij. Experimental decoherence rates are of- 
ten much larger than theoretical estimates com- 
plex systems like conductors or superconductors, many 
low-energy excitations apart from the usual "oscillator 
bath" modes 0, can cause phase decoherence. These 
include charged defects, dislocations, paramagnetic and 
nuclear spins, external flux and voltage noise, and junc- 
tion resistance fluctuations. The advantages of solid- 
state systems (stable circuits, scalability, etc.) cannot 
be exploited until decoherence effects are suppressed. 

It is often assumed that this is purely an engineer- 
ing problem. However this is not true- even completely 
pure systems have intrinsic decoherence, and then theory 
is essential, to see how this can be quantified and con- 
trolled. Here we consider a very interesting possibility, in- 
volving nanomagnetic insulators. These are often chemi- 
cally reproducible, so impurities and defects are kept to a 
minimum- and electronic decoherence is completely ab- 
sent. We show they also have another crucial feature, 
viz., the decoherence can be tuned down to very low val- 
ues. These results are of wide interest, in the search for 
magnetic materials showing spin coherence phenomena. 

1. Intrinsic Decoherence in a nanomagnetic 
" qubit" : At low T the spin Hamiltonian of many large- 
spin nanomagnetic systems (magnetic molecules, rare 
earth ions, or nanomagnetic particles) reduces to a 2- 
state form T-Lo{t) = (Aq-Tj; -|- CqTz), with the Pauli spin 
f acting on the 2 lowest spin levels 0, . The spin gap 
Eq to the next levels is typically 5 — 10 i^, and the 
2-state picture is valid at energies ^ Eq. At higher en- 
ergies the individual electronic spins {s^} are assumed 
to be locked into a nanospin S (so that S — ^j) by 
strong exchange or superexchange interactions. We as- 
sume henceforth an "easy z-axis" nanomagnet; then the 
'bias' energy — gfisSzH^. When €o — 0, the splitting 



Ao between the 2 "qubit" states |— ),|+) (bonding and 
anti-bonding eigenstates of TLo{t)) is produced by tunnel- 
ing between 2 potential wells, each well having a "small 
oscillation" energy Q.o\ typically fio ~ Eq. We define the 
states I t), I i) (eigenstates off,) by |±) = [| T)±| i)]/V2. 
If the total nanomagnetic spin S is not too small, these 
states correspond roughly to semiclassical spin coherent 
states having orientations y\.„ (here <t = which 
depend on both the internal anisotropy field of the nano- 
magnet, and any transverse external field H^. The split- 
ting Ao depends sensitively on H^. 

The intrinsic decoherence in insulating nanomagnets 
comes from entanglement of the nanomagnetic spin wave 
function with that of the nuclear spins and phonons 'lO| . 
Both of these couplings arc well understood. The nuclear 
spins {Ifc} couple to the electronic spins {sj} in S via in- 
dividual hyperfine couplings ^^^s"/^, and the phonons 
couple to S via magnetoacoustic (spin-phonon) interac- 
tions 0. Without yet specifying the precise form of the 
hyperfine couplings (which of course depend on the sys- 
tem being studied) we can define quite generally a vector 
Lj\ which specifies the net effect on the fc-th nuclear spin 
Ife , of all the hyperfine fields coming from the individual 
electronic spins in S: 
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A%{{s'^Y ~ {s'^)') (1) 



Here 4 is a unit vector in the direction of this total field, 
with components , and (s")'^ is the expectation value 
of Sj when S Sn°'. We see that when S fiips from 

Sn-i to Sni, the energy change is just 2wj!4. • Ifc, ie., we 

have a diagonal coupling Tzijj\lk ■ Ifc between the qubit 
and Ifc. One can also write this coupling as f,^,, where 
^2 — (J^lk ■ ifc acts as an extra bias field, in addition to 
the external bias field eo. Note that the weak interaction 
between nuclear spins (typically dipolar) causes slow spin 
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diffusion between them. The net effect of this is to make 
^2 — > ^zit), ie., the nuclear bias field acting on the qubit 
fluctuates in time even when the qubit is frozen. We shall 
argue below that we can neglect this fluctuation in strong 
transverse fields. 

The Ifc also couple to the external field Hq, with Zee- 
man coupling uj^ ihj, ■ Ij. , where 



(2) 



and mfe is a unit vector along Ho. 

If we now take these terms together, we can write the 
interaction Hamiltonian between the nanospin S and the 
nuclear spins {1^} in the form TIms = T^o{t) + V{T,lk), 
where 



V ^ 



k 



Ik 



k 



Wfc ruk ■ Ik 



(3) 



This form is particularly useful for quantifying the deco- 
herence from the nuclear spins. When we come to par- 
ticular examples we will specify the couplings in ©. 
The magnetoacoustic interaction between the qubit co- 
honon coordinate Xn is dominated in 



ordinate f and the 
nanomagnets Ifllllfl 
This term has strength 



by a non-diagonal term Tx J2n 



Cq Xq 



(4) 



where 6'd is the Debye energy, LOq = qcs and Cs is the 
sound velocity. 

Our basic idea is as follows. The allowed nuclear spin 
bath states, when the qubit is in some given state, have a 
density of states which typically has Gaussian lineshape 
[llj l. with a halfwidth Eg; in terms of the {a;|} defined 
above, Eq is given trivially by 



Ei 



(5) 



On the other hand the acoustic phonon energy scale is 
the Debye energy 9d- Now in a nanomagnetic system 
Eq/Od can be < 10^'', suggesting the following tactic for 
suppressing decoherence. If we tune Ao so that 9d ^ 
Aq ^ Eo, then we will be in a "coherence window", 
in which decoherence will be at a minimum because the 
qubit dynamics is too slow to disturb most phonons, but 
too fast for the nuclear spins to react. 

Decoherence Rates: To substantiate this idea, we gen- 
eralise the low field (Aq < Eo) calculations of nanomag- 
netic dynamics |lfll |. where incoherent tunneling relax- 
ation is found, to the high- field regime Aq 3> Eo- Be- 
cause at high field the couplings {^^-"1} of the nuclear spins 
to the qubit are ^ the nuclear Zeeman couplings {w^}. 



and also {w 
turbatively 



Plii 



this dynamics can be solved per- 
We expand about the bare qubit 



Hamiltonian (PJ to 2nd order in (^f/A (assuming eo = 
for simplicity) to get 

Hms = [AoTx + ^ c^fc ?rife • Ik] 



^ ' —{Ik ■ lk){lk' ■ Ifc 



^ 2Ao 

kk' 



(6) 



The decoherence time is defined as the character- 
istic time for decay of the off-diagonal density matrix 
element, starting in state | t)- In the present case we cal- 
culate this matrix element as a path integral over pairs 
of qubit trajectories Tz{t),Tz{t') (each taking values ±1, 
with occasional flips between these values), weighted by 
an influence functional F [Tz(t), Tz(t')] which incorporates 
the interactions 0] . The contribution from the 2nd term 
in ^ to this functional is 0, 0| : 

1^^ = -E^I fdse^-^^[Tz{.s)-4m' (7) 

k 

From this result we can then use standard techniques 
developed for the spin-boson model |^, to find a contri- 
bution 7^ to the dimensionless decoherence rate 7^ = 
I/t^Ao (the commonly used 'decoherence quality factor' 
III is just = 'nh^). We find 



7; 




{ik^\){ik' ^\)^Wlikik' 



9Ikh' 



2A2 



(8) 



which decreases rapidly with increasing qubit operating 
frequency Aq. 

There are 2 other contributions to 7^ coming from the 
nuclear spins First, in writing we omitted a 

renormalisation of the tunneling matrix element caused 
by the coupling to the nuclear spins tl3|, which in fact 
describes the nuclear spin transitions induced directly by 
electronic spin fiips. This adds a contribution 7^ to 7^, 

given by 7^ = 5 Efc |afePi where \dk\ = 7r|wJ!|/2fio. How- 
ever when Ao > Eo, the ratio 7^/7^ 0{Al/nl) < 1 
(the usual WKB reduction of the tunneling amplitude) , 
ie., this term can always be neglected to first approxi- 
mation. Second, we have neglected the instrinsic nuclear 
spin diffusion dynamics, caused by internuclear interac- 
tions. In low fields, when the tunneling is slow, and 
Ao ^ Eo, this intrinsic nuclear dynamics renders the 
tunneling dynamics incoherent [ifl |2j| . However 

when Ao 3> Eo, the nuclear fluctuations are very slow 
compared to the qubit dynamics (of frequency Aq); they 
then add a 'noise' contribution 7^ ^ N/TrAoT2 to 7^, 
where is the number of nuclear spins in each molecule 
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lia. Typically T2" 



10 — 100 Hz at low T in magnetic 

GHz 



molecules 117|), whereas we are interested in Aq 
(see below); thus 7^ will be very small. 

Finally we include the phonon contribution to 7^ . The 
solution to the spin-boson problem for non-diagonal cou- 



pling to phonons gives a contribution 7?'' of form ^ 



7^ = [{S^o^of/^D] COth(Ao/fcBT) 



(9) 



which increases rapidly with Ao. 

The phonon and nuclear spin decoherence mechanisms 
act independently- thus we can get a simple estimate for 
the optimal decoherence rate 7™*" by summing the two 
dominant contributions JHJ and , and minimizing their 
with respect to Ao, assuming ksT < Aq. 



sum + 7^ 
This gives 



ph 



7^^" « V2Si}oEo/9l 
at an optimal tunneling splitting A"^'*: 



(10) 



(11) 



We see that decoherence is minimized for a given S by 
making Eo and flo small, and 0d large, within the con- 
straint that ^ Ao > kBT. If fc^T > Ao we get a 
different (less favorable) answer. 

These simple results actually give reasonably accurate 
results when compared with numerical calculations on 
real systems, as we now see. 

2. Three Examples: We present quantitative results 
for the decoherence rates in 3 different materials. We 
give most details for the Fe-8 example, to illustrate our 
method. 

(i) The Fe-8 molecule: This well-characterized 
molecule Q behaves below ~ 10 if as an electronic 
spin- 10 system, with biaxial effective Hamiltonian 0,0 
HoiS) ^ ~DSl + ESl + Ki{SX + St) - ffMsH^ • S, 
with D/ks = 0.23 K, E/ks = 0.094 K, and K^/kB = 
—3.28 X 10~^ K. The small oscillation frequency fio = 
2SC4VDE - 4.6 where C4 = 1.56 includes the effects 
of the K4 term; T-independent tunneling dynamics ap- 
pears below ~ 0.4 K. The tunneling amplitude Ao(H^) 
is trivially determined by diagonahsation of Ho{S) (Fig. 
1), as are the spin orientations defined previ- 

ously. At a critical transverse field He the barrier is 
destroyed and , n-'' merge. When H;^ is along the 
hard x-axis (so the azimuthal angle ip = or tt), one 
has g^sHc = 2S{D + E) for this Hamiltonian, giving 
He ^ 4.8 T. We find by comparing exact diagonahsation 
with semiclassical calculations, that the latter are accu- 
rate up to ~ 3.3-4.3 T (for = 90°, 0° respectively). 

The hyperfine interactions between the 8 Fe^^ (spin 
5/2) ions and the 205 nuclear spins in the molecule (213 
if ^^Fe isotopes are substituted for ^^Fe nuclei) are of 



2 kinds. The Fe electronic spins interact with any ^^Fe 
ions via contact hyperfine interactions, which we assume 
to be the same as for Fe^^ ions in similar materials [31 . 
On the other hand the hyperfine couplings A^^^s"!^ to 
all the other nuclear spins are thought to be dominated 
by purely dipolar terms [l^ . 




FIG. 1: The tunneling splitting Ao(Hx) for the Fe-8 
molecule, for transverse field Hx in the xy plane, at azimuthal 
angle (j> from the in-plane hard axis. We use parameters given 
in the text 0,0 for the molecular spin Hamiltonian 7io(S). 

Using the known nuclear positions and moments, the 
hyperfine interactions and the nuclear Zeeman couplings 
are then quantified numerically, to find the {^^-"1} and 
hence Eg (see Fig. 2), and also the {w^}. In spite of 
the large number of protons, Eg is quite small (partic- 
ularly when we substitute ^H nuclei for ^H), making 
_Fe-8 a reasonable candidate for coherent dynamics in 
high fields. One may also calculate the nuclear dynam- 
ics 0], but here we simply note that NMR experiments 
[13 indicate that at low T, T2 ~ 10 - 30 msec for ^H 
nuclei in magnetic molecules. We then calculate numer- 
ically the different nuclear spin contributions to 7,^ (ie., 
the dominant contribution JH)) , plus the smaller contribu- 
tions 7^ and 7^^), and also the phonon contribution 7^'', 
as a function of H;^. 

The results are shown for high fields in Fig. 3. Actually 
almost all experiments on the quantum dynamics of Fe- 
8 have been done in the regime Aq < i?o, ie., where we 
expect nuclear spins to cause incoherent tunneling llOl- 
which is what is found experimentally 17., 2M, ■ 
However once Aq exceeds Eo, all nuclear spin decoher- 
ence should fall off very fast. The nuclear noise contribu- 



tion 



< 10-^ once |H^| > 3.2 T (for = 0°) or > 2 T 
(for (j> > 90°), so we can safely ignore it. The estimate 
(Uni then gives 7^™ « 5.8 x 10"^ at a A°p* « 0.14 K, 
when (p = and H^ = 3.45 T (assuming the optimal set 
of nuclear isotopes, with ^H instead of ^H, no ^^Fe, etc). 
In Fig. 3 we calculate 7^ numerically, adding all nuclear 
spin and phonon contributions. For = we get a nu- 
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merical value 7^™ = 6 x 10-^ at a A°p* = 0.135 K when 
= 3.45 T; thus the estimate works well for (f) — Q. 
The further reduction of 7™™ to ^ 2.4 x 10^^, when 
(j) 90" {H^ — 2.93 T), conies from decreases in both 
flo and the effective S (since both the tunneling barrier 
and the change jS^ — decrease). Fig. 3 clearly illus- 
trates the window for coherent dynamics which opens up 
at high transverse fields {H^ ^ 2.9 — 3.4 T, depending 
on <b) in Fe-8. 
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FIG. 2: A histogram of the calculated couplings {ujI} be- 
tween the -Fe-8 spin-10 and both ^ H nuclei and ^"^ Fe nuclei, 
binned in 0.1 MHz intervals (the couplings to the other nu- 
clear species are not displayed). Apart from the core polar- 
isation coupling to ^""^ Fe nuclei, all hyperfine couplings are 
assumed dipolar, and calculated using the known positions of 
each nucleus [231. We assumed Hartree-Fock wave- functions 
for each Fe'^^ ion (correcting the results of ref. |25ll . which 
used an incorrect coordinate system for the molecule). The 
insert shows the variation with of the nuclear multiplet 
linewidth parameter Eo (see text). Curve (2) has naturally 
occurring isotopic concentrations. Curve (1) substitutes ^"^ Fe 
for ^'^Fe, and curve (3) substitutes for ^H. Previous mea- 
surements of the "hole- width" parameter ^o, which has Eo as 
an upper bound (taken from ref. [T5|| ') are also shown for 2 of 
these cases. 



(ii) A^i-based particles: Consider now a pure Ni 
particle at low T, with Qq = 1 K and spin S, coupled 
to a substrate with 6'£) = 300 K. The concentration 
of ^^Ni nuclei (spin- 1/2) is 2:59 ^ 0.01 in natural Ni, 
and the hyperfine coupling uj. Uo ^ 1.4 mK. From 
(Uni we have 7^" - SEo/lO'^ - OiS^/^tOoxH^ /lO^), 
since Eo ^ ujoN^^'^, where ~ S" is the number of nu- 
clear spins. Then in natural Ni, 7™™ reaches unity once 

5 - 2.6 X 10-*, for a A^"^*^ - 0.2 K. Isotopic purifica- 
tion by a factor 10^ (a major undertaking!) would reduce 
7™™ by only a factor of 10. Note that Ni has an unusu- 
ally low concentration of nuclear spins, with rather weak 
hyperfine couplings. Thus this example teaches us that to 
see any macroscopic coherence effects in tunneling mag- 
nets for S > 10^ will require almost complete isotopic 



purification (there are no known magnetic species not 
having at least one natural isotope with non-zero nuclear 
spin). Note further that we only consider here intrinsic 
decoherence from phonons and nuclear spins inside the 
particle- we have ignored electronic decoherence, which 
certainly exists in Ni. 




Hi(T) 



FIG. 3: Calculated behaviour of 70(Hx) for an Fe-8 molecule, 
using the results for Ao{ili', the {ljJJ(Hx)}, and Eo (Fig. 1), 
and also evaluating f2o(Hx) as a function of Hx. We show the 
individual contributions from 7^ for different field orientations 
in large transverse fields, for 7^ (shown for cj> — O"), and 
for the phonon contribution 7^'', assuming a Debye energy 

Od — 33 K (the small contribution 7^ is not shown). Results 
are shown for the optimal distribution of nuclear isotopes, ie., 
with '^H substituted for ^H, and only ^Ve, ''^Br, "iV and 
^^O species in the molecule. 



(iii) Ho ions: Not all systems will have an opti- 
mal intrinsic decoherence described by (|10(l . Consider 
the LiHoxYi-xF4 svstem |2(Tj |. in which sharp absorp- 



tion lines are seen 



system j. 



at low X. The Ho^ ions 



(with spin = 8) have a lowest doublet state, with split- 
ting Ao oc |H^p in a transverse field- only one other 
state, at an energy f^o ~ 10.6 K, is important for the 
low-energy physics (all other electronic spins levels are 
at energies ^ 100 K). The new feature here, which ren- 
ders H10|l inapplicable, is a very strong hyperfine coupling 
0.039 K to the / — 7/2 Ho nuclear spin- whereas 



-'Ho 



the couplings to the "satellite" Li, Y , and F nuclear spins 
are very weak. Thus instead of having a Gaussian line- 
shape, the hyperfine multiplet has a "toothcomb" struc- 
ture with 15 Ho lines spaced in intervals of 0.039 K , 
each weakly broadened (by less than 1 mK) by the other 
nuclei. The quantity Eo is then not well-defined [ll| . 
However one may instead just calculate the coupled Ho- 
nuclear spin dynamics in 2nd-order perturbation theory, 
since ilo 3> all hyperfine couplings, and all Ho transi- 
tions go via the single intermediate level. We then find 
the optimal strategy is to (a) freeze the Ho nuclei by cool- 
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ing to ksT <C 0.039 K, and (b) make H;^ large enough 
so that precessional decoherence from the satehite nuclei 
is eliminated, but small enough so that Ho nuclei are 
weakly excited. We find that at T = 3 mK, tuning of 
Ao to 30 mK (using ~ 0.3 T) gives a contribution 
7™'" ~ 1.8 X 10~^ from the Ho nuclei (phonons giving 
a contribution < 10~^ here). Some residual decoherence 
also comes from the 'satellite' nuclear spins- in this con- 
text it is interesting that Ho~^^ ions can be prepared in 
hosts with almost no other nuclear spins (eg., CaWO^, 
with isotopically purified W\ see [13 )■ Thus we can ex- 
pect a reduction of 70*™ to very low levels, limited only 
by a weak phonon effect, external noise, and any spin 
impurities. 

3. Discussion: In insulating systems where nuclear 
spin decoherence is suppressed by a transverse field, de- 
coherence optimisation is well described by (|10() . This 
means having a small S, small i?o, and a "stiff' system 
(high 9d)- Surprisingly it also means low r^o, ie., weak 
magnetic anisotropy. By taking all these measures, very 
low values of 7^ can be attained. 

Our results are thus good news for coherence in small 
nanomagnets, and magnetic qubits. In most materials 
a correct choice of parameters requires strong transverse 
fields, but in systems like LiHoxYi-xFi with one very 
strong hyperfine coupling wjj , and otherwise weak or non- 
existent ones, it pays to keep Aq low (ie., not too strong 
fields), and to have kT <C (to suppress thermal nu- 
clear spin noise). 

For genuinely macroscopic superpositions of magnetic 
states, the result ()10(l is not such good news. Although 
one may stave off decoherence for large S by isotopic 
purification of nuclear spins, our Ni example shows this 
tactic reaches its limit once S ^ 10^ — 10^. 

So far experiments on nanomagnets have concentrated 
on field ranges where Aq < i?o, and so incoherent tun- 
neling is observed- any experiments on, eg., in the 
range |H^| '--^ 3 — 4 T would be of great interest. In al- 
most all experiments so far, on molecules or rare earth 
ions, collective tunneling caused by inter-spin dipolar in- 
teractions complicate the interpretation 20, 21, 22] . A 
clean observation of coherence, with a measurement of 
70, will thus involve manipulation at microwave frequen- 
cies, on a properly isolated single system. 

We thank G.A. Sawatzky for useful discussions, and 
NSERC, the CIAR, and grant number NS-1767.2003.2 
in Russia, for support. 
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